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On ' m „ 

In the heavy quark limit, the heavy baryons Uq (lo stands for £, S or SI 

Ph. and Q = b or c) are regarded as composed of a heavy quark and an axial 

vector, light diquark with good spin and isospin quantum numbers. Based 

on this diquark picture we establish the Bethe-Salpeter (B-S) equation for 
• , , 

lOq in the limit where the heavy quark has infinite mass, tuq — > oo. It is 

found that in this limit there are three components in the B-S wave function 

(*) 

for 0Jq . Assuming the kernel to consist of a scalar confinement term and 
a one-gluon-exchange term we derive three coupled integral equations for 
the three B-S scalar functions in the covariant instantaneous approximation. 
Numerical solutions for the three B-S scalar functions are presented, includ- 
ing a discussion of their dependence on the various input parameters. These 

solutions are applied to calculate the Isgur-Wise functions £(u>) and £(u>) 

(*) (*) 

for the weak transitions ~~ > • Using these we give predictions for 
the Cabibbo-allowed nonleptonic decay widths and up-down asymmetries 
for $7;, — > plus a pseudoscalar or vector meson. 



PACS Numbers: ll.10.St, 12.39.Hg, 14.20.Mr, 14.20.Lq 



I. Introduction 



The physics of heavy hadrons has been a subject of intense interest in recent 
years. One reason for this is that more and more experimental data are being 
accumulated. Another reason is the discovery of the new flavor and spin symme- 
try in QCD, SU(2)f x SU(2) S , in the heavy quark limit and the establishment of 
heavy quark effective theory (HQET) PL However, in comparison with the heavy 
meson case, heavy baryons have been studied much less, both experimentally and 
theoretically. 

On the other hand, the experimental situation with heavy baryons has been im- 
proving recently, with more measurements becoming available. For instance, OPAL 
has measured some physical quantities for A b , such as its lifetime and the production 
branching ratio for the inclusive semileptonic decay A b — > Al^uX ||. Furthermore, 
measurements of the nonleptonic decay of A& have also been made, through the well- 
known process A& — > AJ/ip. The discrepancy between the measurements made by 
UA1 || and those by CDF H and LEP || appears to have been settled by the new 
measurement from CDF M. However, compared with D and B mesons, the data for 
heavy baryons are still very limited. Besides A b , there has been few data on other 
bottom baryons |/J, although we expect more data to appear in the near future. 
Clearly the time is right for serious theoretical studies of heavy baryon properties 
to begin. 

Theoretically, the HQET can simplify the physical processes involving heavy 
quarks, since with the aid of the HQET the number of independent form factors is 
reduced. For instance, in leading order of the 1 / tuq expansion only one form factor 
(the Isgur-Wise function) remains for the A& — > A c transition, while for — > u^ 1 

(we follow the notations of Ref. ||, u could be S, S or Q and Q = b or c) there are 

(*) 

two independent Isgur-Wise functions. (Note that ujq is a notation implying either 
ujq or Uq.) The behaviour of these functions depends on the nonperturbative effects 
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of QCD which control the dynamics inside a heavy hadron. Hence some nonpertur- 
bative QCD model has to be adopted from which these Isgur-Wise functions can be 
obtained. In previous work ||, we established the B-S equation for Aq, which is as- 
sumed to be composed of a heavy quark, Q, and a scalar diquark. Some theoretical 
predictions for A& — ► A c were also obtained. It is the purpose of the present paper 
to generalize such an approach to the heavy baryons, 0Jq , and consequently give 
some phenomenological predictions for the weak decays of such baryons. 

When the quark mass is very heavy compared with the QCD scale, Aqcd, the 
light degrees of freedom in a heavy baryon Aq (Q = b or c) become blind to the 
flavor and spin quantum numbers of the heavy quark because of the SU(2) f x SU(2) S 
symmetry. Therefore, the light degrees of freedom have good quantum numbers, 
including angular momentum and isospin. These quantum numbers can be used to 
classify heavy baryons. For example, the light degrees of freedom of Aq have zero 
angular momentum and isospin. For Eq\ the angular momentum and parity J p of 
the light degrees of freedom are 1 + , and the isospin is also 1 in order to guarantee 
that the total wave function of the hadron is antisymmetric. Hence it is natural to 
consider the heavy baryon to be composed of a heavy quark and a light diquark. 
This is our underlying assumption. 

Based on the picture of the composition of the heavy baryon which we have just 
presented, the three body system is simplified to a two body system. We establish 
the B-S equation for the heavy baryons, Uq , in this picture. The heavy quark 
symmetry can be used to simplify the form of the B-S wave function greatly. It can 
be shown that in the limit rriQ — > oo there are three components in the B-S wave 
function, and hence we have three corresponding scalar functions. We solve the B-S 
equation numerically by assuming that the kernel contains a scalar confinement term 
and a one-gluon-exchange term. The explicit dependence of the B-S wave function on 
the parameters of the model will be discussed. Furthermore, we calculate the Isgur- 
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Wise functions in terms of the B-S wave functions and give theoretical predictions 
for the Cabbibo-allowed two body nonleptonic decays Clb — ► il^ plus a pseudoscalar 
or vector meson. 

The light degrees of freedom of uj^q belong to a 6 representation of flavor SU(3). 
Taking Q = b as an example, a>j includes Ej,*^' ' - , and fl^ ~. The total 

spin of ujq and Uq are | and | respectively. There is no strange quark in t!q while 
there is one strange quark in Eg and two in VLq respectively. 

The remainder of this paper is organized as follows. In Section II we establish the 
B-S equation for the heavy quark and axial vector light diquark system and discuss 
the form of the kernel. In Section III we derive explicitly the coupled integral equa- 
tions for the B-S scalar wave functions. In Section IV we discuss the normalization 
condition of the B-S wave function by exploiting the normalization of the Isgur-Wise 
function at the zero recoil point. The numerical solutions of the B-S equation and 
their dependence on the parameters in our model are presented in Section V. In 
Section VI we calculate the Isgur-Wise functions and give predictions for the decay 
widths and up-down asymmetry parameters for flj, — > f)W plus a pseudoscalar or 
vector meson. Finally, Section VI contains a summary and discussions. 

(*) 

II. The B-S equation for luq 

As discussed in Section I, uJq 1 is regarded as a bound state of a heavy quark, 
ipQ, and a light axial vector diquark A^. In the following, uq denotes the Dirac 
spinor of ipQ and 77^ represents the polarization vector of A^. Then = uqT]^ 
can be decomposed into spin-| and spin-| states [TI], |TT| which represent ujq and 



ujq respectively. These two states are degenerate in the heavy quark limit. Using 
the notation of Refs. [|TT| [12| this doublet is described by B™(v), where m — 1,2 
correspond to ujq and Uq respectively, is the velocity of the heavy baryon and 



3 



Bf, = B^ + [0. Explicitly, we can write 

B U V ) = ^(7m + v a«)T5«(v), 

Bl{v) = «», (1) 

where u(v) is the Dirac spinor and u^{v) is the Rarita-Schwinger vector spinor. 
B™{y) satisfies the following conditions: 

fB?(y) = B?(v), v*B?(v) = 0, = 0. (2) 

The above constraints for m = 1 can be seen from fu{v) = u(v) while for m = 2, 
they are the properties of a spin-|, Rarita-Schwinger vector spinor. 



Under a Lorentz transformation A 1C 



B, -> A^D(A)B„, (3) 

where -D(A) is the spinorial representation of A. Under the heavy quark spin 
transformation II 

£ M - -T5#^, (4) 

where e = ei,ez,es are three mutually orthogonal four-dimensional unit vectors 
which are also orthogonal to v, (i.e., e* • v — 0), and = — 1 (z = 1, 2, 3). The three 
unit vectors are associated with the heavy quark spin operators which are generators 
of the SU(2) S symmetry. 

Since uuq is composed of ipg and A M , we can define the B-S wave function of 

(*) u 
u q by 

Xfl (xi,x 2 ,P) =< 0\T4j q (x 1 )A i1 (x 2 )\J* ) (P) >, (5) 

where P = m^(*)V is the total momentum of Uq and v is its velocity. Let mq and 
mo he the masses of the heavy quark and the light diquark in the baryon. Let us 
define Ai = — ^ — and A2 = — ^P- — and let p be the relative momentum of the 

1 mQ+m D A m Q +m D r 



two constituents. The B-S wave function in momentum space is denned as 

X>i,x 2 ,P) = e iP X J 7^^- X ™>), (6) 

where X = XiXi + X2X2 is the coordinate of the center of mass and x = x\ — X2- 
The momentum of the heavy quark is pi = X\P + p and that of the diquark is 
p 2 = -A 2 -P + p. 



It can be shown that Xp m {p) satisfies the following B-S equation fT3[ 

X Pm {p) = S F (X 1 P+p) J ^G pu (P,p,q)x P Jq)S%>(-\ 2 P+p), (7) 



where G pu (P,p,q) is the kernel, which is defined as the sum of all the two particle 
irreducible diagrams with respect to the heavy quark and the light diquark. For 
convenience, in the following we use the variables 

pi=v-p - A 2 rry*) , p t = p - {v ■ p)v. (8) 

Q 

Then in the leading order of the 1/rriQ expansion we have 

S F (X 1 P + p) = l ^ + J ] (9) 

2(p l + E + m D + le) 

where E is the binding energy. From Eqs.(^) and (Q) it follows that Xp m (p) satisfies 
the following equation 

^XPm(p) = XPrniv)- (10) 

The propagator of the light axial vector diquark has the form 

S D „„te) = - V-^f>° (ll) 
P2 — m z D + le 

In the limit tjiq — > 00 we have p 2 = —rri£)V + p and hence we have 

S»M - "» pf-W+ie ' (12) 



where W p = yp* + m D- The corrections to Eqs.(^) and ( jl^D are from 0(l/mg) 
terms. 

Now we discuss the form of the B-S wave function xp m {p)- I n the heavy quark 
limit, due to the 577(2) s x SU(2)f symmetry, the internal dynamics of the heavy 
baryon, ujq , is determined by the light degrees of freedom and the flavor and spin 



direction of the heavy quark, Q, is irrelevant. Consequently we havefllQ, 12 



Xp(p)=u Q (v)r)„(^(v,p), (13) 

and 

<0|A M |light,l + >=^C M >,p). (14) 



Since v ■ rj = 0[|10), the tensor (^ u (v,p) can be expanded as 

C»(v, p) = A l9 ^ + A 2 v»p v + A^p v , (15) 

where Ai(i = 1,2,3) are Lorentz scalar functions. After expressing UqT]^ in Eq.(O) 
in terms of B™(v) (m = 1,2), we have the following form for the B-S wave function 

XPm = A x B^{v) + A^p v B v Jv) + A^p v B v Jv). (16) 

Therefore, we have three components in the B-S wave function, xpm{p)> an d they 
correspond to three scalar B-S functions A{(i = 1,2,3). This is consistent with our 
diquark picture for lOq 3 . In the heavy quark limit, the dynamics inside the heavy 
baryon is controlled by the configuration of the light degrees of freedom. Since the 
light diquark is a 1 + object, it has three different configurations. Consequently there 
are three components in the B-S wave function which describe the dynamics in the 
heavy baryon ujq . 

In fact, we can derive the form of Xp m {p) i n another way. We may first write 
out all the possible terms which have the same behaviour as Xp m {p) under Lorentz 
transformations. Then by applying the condition Eq.fllOD and ensuring the proper 



behaviour under the heavy quark spin transformation, Eq.(|j), we obtain the same 
result as given in Eq. (|T6D . 

Considering p M = v ■ pv^ 1 + pf, and using the constraint v^B™^) = 0, it will be 
convenient to define 

A = A U C = A 2 + v-pA 3 , D = A 3 , 
which results in the following expression for the B-S wave function: 

X £ m = AB^v) + Cv» Vtv B v m (v) + Df tVtv B v m {v). (17) 

A, C and D in Eq.(|P7|) are functions of pi and p\. Their behaviour is controlled 
by nonperturbative QCD. Our aim is to obtain explicit forms for them with some 
QCD-motivated model for the form of the B-S kernel. 

Motivated by the success of the potential model |14j], scalar confinement and one- 



gluon-exchange terms were used in the kernel when studying Aq in Ref. 0. This 
form was also used in the heavy meson case in Ref . [T5| . In the present work we will 
also adopt this form of the kernel 

iG pu = gP uj g IVi +VfM0 T^ v V 2l (18) 

where T^ pu is the vertex of a gluon with two axial vector diquarks. This vertex 
should reflect the internal structure of the diquark. In this work, we use the model 
proposed in Ref. ]TB| where this vertex has the following form (see Fig. 1) 

- l ^ 9s Y^F v {Q 2 ) ) (19) 

with 

iw = ( P2 + v ' 2 y g v e - (p&r + pgy 1 " )■ 

In Eq.([l9|) g s is the strong interaction coupling constant and Fv(Q 2 ) is introduced to 
describe the internal structure of the axial vector diquark. The form factor, Fy(Q 2 ), 



depends on nonperturbative QCD interactions and will be determined phenomeno- 
logically, by comparison with experiment. 

As discussed in Ref. [f|, when we consider the vertex of two heavy quarks with 
a gluon, the momenta of the two heavy quarks are p\ = Xim^ Q v + p and p[ = 
^i m AQ v + Q respectively, where p and q are relative momenta and of the order 
Aqcd- in the heavy quark limit the heavy quark is almost on-shell and moves with 
constant velocity. It can be shown that p\ = q\ at this vertex when the heavy quark is 
exactly on-shell. This is the so-called covariant instantaneous approximation |], [15]. 
With this approximation, V\ and V 2 in G pu (P,p,q) are replaced by 

7 < = V i | B=aft (i = l,2). (20) 



III. Coupled integral equations for three B-S scalar wave 
functions 

In this section we will derive explicitly three coupled integral equations for the 
B-S scalar wave functions. Substituting Eqs.@ and (|l|) into Eq.(|7|) and considering 
the form of the kernel in Eq.(|l8"D and the property in Eq. fllCf) , we obtain the following 
form for the B-S equation 

[pi + E + m D + it) [pf - + le) 



where 

(v^Pp + p^Vp) p^p p 



M mfl = i 



9^ P - v^Vp + 



m D m D 



iV<l [G^(P,p,q)xPmM} 



(2vr 



Pi=ni 
(22) 



and we have made explicit use of the covariant instantaneous approximation. 

Substituting Eqs.(|l7D, flUf ) and (|19D into Eq . (p2|) and again using the covariant 
instantaneous approximation we have 

= B^v) J ^A(V X + 2 Pl V 2 ) + ±-vPj ^{ Pt ■ B m (v)A(Vi + 2 Pl V 2 ) 
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Hvi + m D )[-Apt ■ B m (v)V 2 + q t ■ B m (v)(CVi - Dp t ■ q t V 2 )] 

d 4 g 

{ Pt ■ B m {v)A{V x + 2 Pl V 2 ) +pi[-A Pt ■ B m (v)V 2 



+p t ■ q t q t ■ B m (v)[-CV 2 + D(V X + 2 Pl V 2 )]} - \p» f 

771 j-) J 



(2tt)< 



+q t ■ B m (v)(CVi - Dp t ■ q t V 2 )] +p t ■ q t q t ■ B m (v)[-CV 2 + D{V X + 2 Pl V 2 )]} 
+ / -Mid* ■ B m {v)[-CV 2 + D{V X + 2 Pl V 2 )). (23) 

We notice that in Eq. (p3|) there are terms of the form J j^Ssqtf and / -SS^q^tf^ 
where / is some function of p 2 , q 2 , and p ■ q. On the grounds of Lorentz invariance, 
in general we have 

<V<I tff = fiv" + f 2 pt, (24) 



(27T) 



and 

d 4 g 



/ j2^ q ' qU = 9xfV + 92V " VU + 9 ^ ft + 94VUp " + 95P " P "- (25) 



From Eq.(J2J) only the f 2 term can contribute, while from Eq.(Pq) only the gi,gs 



and g§ terms can contribute, since v v B v m [v) = 0. It can be easily shown that 

r d 4 q pt-qt . 
12 J (2vr) 4 Vl ^ 

f d 4 g (pt ■ q t ) 2 - p 2 t q 2 t 
91 J (2vr) 4 2p 2 L 

93 = 0, 

f d 4 q 3(p t ■ q t ) 2 - p 2 q 2 
95 ~ Jj2^ 2rf ^ (26) 

With the aid of Eqs. fl2"4|) , (25|) and (|26|) we can express M£ in Eq . (|23"D in terms 
of B^(v), v^ptvB^Jyv) and p^pt v B v m [v) . Let us define 

A(P 2 ) = J ^(Pi,P 2 t ), C(p 2 ) = J ^ C Ml Dip 2 ) = J ^D( PhP 2 ), 

(27) 

where A, C and D are functions of p\ only. Then one obtains the expression 
M£ = B*(v) J || Ufa + 2 Pl V 2 ) -C ^- ~ P ^ V 2 



+o ( *-*'V'ft ? (Vi + W)" 

^Pt 

1 ,. „ , , f d 3 q t 



m 



-v^pt ■ B r , 



v 



-C 



(Pi - m D ) 



2 \ Pt ft T ~ 



Pt 



Pt 



Vo 



+D 



nipt ■ Qtf 



Pt 



Wi + (pf + m 2 D )V 2 ] 



rn 
+C 

+D 

+Pi 



-pfp t ■ B m (v) 



D 



d?q t 



{AiK + pM) 



PfQuy . m l(3(p t ■ q t ) 2 - p 2 t q 2 t ) + 2p 2 (p t ■ q t f ~ 

Pi —Vl H — V 2 



Pt 2Pt 
m 2 D (3(p t ■ q t ) 2 - p 2 q 2 ) + 2p 2 (p t ■ q t f 



2pt 



(Vi + 2 Pl V 2 ) 



(pt ■ qtf,~ 



pi 



Vo 



(2? 



In Eg . there are poles in pi at W p — ie,—W p + it and — E — mo — ie. By 
choosing the appropriate contour, we integrate over pi on both sides of Eq. (|21~D and 
obtain the following, three coupled integral equations for A,C and D 

-1 r d 3 q t 



MPt) 



2W P (E + m D - W p ) J (2tt 



{A(q 2 )(V 1 -2W p V 2 ) 



-C(q{) 



.Apt-qt?-- 2 " 2 



PtVtTl 



2p 2 t 



V2 + D(qf) 



(pt -qt) 2 - p 2 q 2 

2p 2 t 



{Vi - 2W P V 2 ) 



(29) 



C(P 2 ) 



d 3 q t 

2m 2 D W p (E + m D -W p )J (2tt) 3 
i(E + m D )W p + m 2 D )V 2 \ 
Pt ■ qt 



{A{q 2 )[W p V x 



+C(qf) 



+Hq 2 t) 



pi 
(pt ■ qtf 
p 2 



((E + m D )W p -mj D )V 1 + W 1 



(pt-qt) 2 ,~ 



V 2 
Pt 



V, 



[-W& + ((E + m D )W p + m 2 D )V 2 ] } , 



(30) 



Dip!) 



d 3 q t 

2m 2 D W p (E + m D - W p ) J (2tt) 3 



{AiqDiK-WpV,) 



-C(q 2 ) 



Pt ■ qt 
Pt 



W p Vi + 



m 2 D (3(p t ■ q t ) 2 - p 2 q 2 ) + 2p 2 (p t ■ q t ) 2 ~ 



2 P f 



V, 
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-D{q 2 t ) 

(Pt ■ qt) 
Pt 



m 2 D (3(Pt ■ qt) 2 ~ P 2 q 2 ) + 2p 2 (p t • qtf 



2pf 



(Vi - 2W P V 2 ) 



-w p v 2 



(31) 



If one knows the form for the kernel, V\ and V 2 , then A(p 2 ),C(p 2 ) and D(p 2 ) can 
be obtained from Eqs.(|29|), ( pCf ) and (|3T|) . Consequently from Eqs. (JT7|) , (f21|) and 
(|28| ) we find the following expressions for A(pi,p 2 ), C(pi,p 2 ) and D(pi,p 2 ): 

d 3 q t 



MPhPt) 



(pi + E + m D + ie) (pj - W 2 + ie) J (2vr 



{i(g?)(Vi+2p,V 2 ) 



-C{q 2 ) 



(Pt ■ qt? - p 2 qU~ 



2p 2 t 



V 2 + D(qf) 



{Pt ■ qt) 2 ~ P 2 t q 2 n ~ 



2p 2 t 



{V x + 2 Pl V 2 ) 



(32) 



C( Pl ,p 2 t ) 



d 3 g t 

m 2 D (pi + E + m D + ie)(pf - W 2 + ie) J (2tt) 3 



{-i(g t 2 )[^ 



+(p 2 + m 2 D )V 2 }-C(q 2 ) 



, > 2 \Pt ■ qt-rr 
(Pi -%) — h—Vt +Pi 



ipt ■ qt) 2 ,~ 



p 2 



v, 



pi 



+D(q 2 ) 



(pt ■ qtf 
Pt 



[ Pl V 1 + (p 2 + m 2 D )V 2 } 



(33) 



D(pi,Pt) 



m 2 D (pi + E + m D + ie)(pf - W 2 + ie) 



^{MtiW+pM) 



+C(q 2 ) 



+D(q 2 



Pt ■ qt T > , m 2 D (3(p t ■ q t f ~ P 2 q 2 ) + 2p 2 (p t ■ q t f ~ 
— —Pi v i H —a V 2 



Pt %Pt 
m 2 D (3(p t -qt) 2 - PhD + 2 Pt(Pt ■ qt) 2 



2pf 



{Vi + 2 Pl V 2 ) 



{Pt ■ qtf 
Pt 



PiV 2 



(34) 



A model kernel, specified in terms of V\ and V 2 , for the B-S equation in the scalar 
light diquark case was given in Ref.|J. In the present axial vector diquark case, a 
model vertex for a gluon with two 1 + diquarks is given in Eq. (|19|) , where Fy(Q 2 ) 
describes the internal structure of the light diquark. Following Ref.[16| we take the 
form for Fy(Q 2 ) as 



F V (Q 2 



a 



Q 2 + Qi 



2 ■ 



(35) 
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where Q\ is a parameter which freezes Fv(Q 2 ) when Q 2 is very small. In the 
high energy region the form factor is proportional to 1/Q 2 , which is consistent 
with perturbative QCD calculations [[HJ. By analyzing the electromagnetic form 
factor for the proton it was found that selecting Q\ = 3.2GeV 2 can lead to results 



consistent with the experimental data |16[. Note that in Eq . (|35|) we do not consider 



the difference between longitudinal and transverse polarization states. The reason 
is that we are considering the bound state with a binding energy of order Aqcd, so 



this difference should be small (see Ref . |16| , where this difference is a factor q 2 +q-2 



with Q\ approximately 15GeV 2 , so this factor is close to 1 in our discussions). 

Based on the above discussion, the kernel for the B-S equation in the baryon 
case is taken to have the following form 

~ 8lTK 33 f d 3 k 87TK 

V i = 77 75 I 9i9 ~ ( 27r ) 6 (Pt ~ Qt) 



i(Pt-qt) 2 + fi 2 ] 2 V ' V " ^ J (2^)3(^ + ^)2' 
16tt a (eff)2 2 

_ 1U7T U. s Vi , , 

2 3 [(pt-qt) a + lf][(pt-<k) a + Cfi\' 1 ) 

where k and af s > are coupling parameters related to scalar confinement and the one- 
gluon-exchange diagram respectively. The second term in Eq. fl36|) is the counter term 
which removes the infra-red divergence arising from the linear confinement in the 
integral equation. The parameter [i is introduced to avoid the infra-red divergence 
in numerical calculations. The limit \i — > is taken in the end. Besides Qf, there 
are two parameters k, and c4 eff \ in the kernel. However, they should be related to 



each other when we solve the coupled integral equations (P9f), (pOf ) and (|31|) . This 
will be discussed in detail in Section V. 

Since we now have an explicit form for V\ and V2 in Eq. (|36D, we can reduce 
Eqs. (p9|) , (p0|) and (0) to one dimensional integral equations. With the aid of the 
formulas given in Appendix A we obtain the following equations from Eqs.(p9|), ( p0| ) 
and (1311). 



A (Pt) = ot jt" rET—— 777T / -^{[^ KF i{\Pt\Mt\) ' 



2W P (E + m D - W p ) J 4vr 2 u 3(Q 2 - /i 2 ) 

12 



(F 2 (\p t \, \q t \,Li) - F 2 (|p t |, \q t \, Q,))]A{q 2 t ) - tornF^], \q t \)A(p 2 t )} 
1 f q 2 dq t 8tt/3 2 

+i ? 2(b*Ug t |,gi)]-4[( F 4(b*l,|g*l,/i)-i 71 4(b*|,|g*|,Qi)]}C'(g, 2 ) 
p< 

+ 2^( g „ + 1 mg -^) /l^ {4TKb?F ' (lKl - W) " 19,1)1 

\Qixf3W p 2 , 

5vWt( F 2(|Pt|, |gt|,/x) - ^(M, M,<2i)) 

"'Wi — ^ ) 

+ ^(F 4 (\p t \,\q t \^)-F 4 (\p t \,\q t \,Q 1 ))]}D(q 2 ), (37) 
^Pt 



(fiflftl, IftU) - fiflftl, l<fcl,<W)((£o + >"d)Hv + m 2 D )]A(<i) 

s™w M M) Am - 2mW ' +MJ _ ff>) /f 

{[8^,(| ftl , l9 ,|) m °- (g ° 2 +T "° ) ^ 

Pi 

1 ^dftrro_«,^i 



16tt/? 



3(Q?-/i 2 



■(■PkdPtl, |&|,A*) - #i(bt|, 



(Eo + rnjWr + ml ^ _ g^^^ | ft | )5(p?)}> (38 ) 

Pi 

(F 2 (\Pt\, |gtU) - ^(bt|, |g t |, Qi))]i(g 2 ) - 87r«Fi(M, M)i(p 2 )} 

+ 9 2H/ , p j ==t / 4%{[-8^F 3 (b t |,|g t |)^ 

2m 2 D W p (E + m D - W p ) J 4tt 2 p£ 

+^?^(^(|Pti, |%U) - F 4 (\ Pt \, \ qt i Ql ))M + *™z 

o\t^ 1 — fi ) zp t 

+8nKW p F 1 (\p t \,\q t \)C(p 2 t )} 
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{[87rKF B (|ptUgt|)- 



q 2 t m 2 D 16vr/3 



W p (p? + 3mf,) 



Pt 



3(Q 2 -/i 2 ) v ^'"^"^ p 2 
-8ttk( P 2 + m^Fid^l, |%|)£(p 2 )}, (39) 



where |p t | = yp 2 . The functions Fj(z = 1,...,5), appearing in Eqs.(|37D, (|38"D and 



(|39|) , are defined in Appendix A. From these three coupled integral equations we 
can solve numerically for A(p 2 ), C(pf) and D(p 2 ). This will be done in Section V. 
In the next section, we will first discuss the normalization of the B-S wave function. 

IV. Normalization for the B-S wave function 

It can be seen that the overall normalization of A(p 2 ), C(p 2 ) and D(p 2 ) cannot 
be determined from Eqs.(^), (^) and ([39]). With the help of heavy quark symme- 
try, the normalization constant can be obtained from the fact that the Isgur-Wise 
function is normalized to one at the zero-recoil point. In the limit mq — > oo the 
weak transition matrix element induced by the current cTb for uj^ — > lu^ has the 
following form from the HQET, 

U*\v')\cTb\J b *\v)) = B^{v')TB^{v){i{u)g, u + C(^H<), (40) 

where to = v • v' is the velocity transfer, m, m' could be 1 or 2, and T is an arbitrary 
Dirac matrix. At the zero-recoil point, v = v', only the ^(o;)^ term contributes 
and we must have 

= 1) = 1. (41) 
On the other hand, the transition matrix element for U: ] - w w is related to 
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the B-S wave functions of uo^ and by the following equation 

(ul*\v>)\cTbU*\v)} = j ^r P Up')rxpJp)s»lM, (42) 

where P (P 1 ) is the momentum of iJf* {^i*^) and Xp'm'ip') * s the wave function of 
the final state u>^'(v'), which satisfies the constraint 

Xp>m>(p')f = X P >m>(p')- (43) 

At the zero- recoil point, p' = p, since the light diquark sees no change in the heavy 
quark part it does not change its relative momentum. 

The scalar B-S functions of the final state B-S wave function obey the same B-S 
equation as fl29|), © and ([31]). Substituting Eq.(0) into Eq.(gD and using Eq.([p 
we have 

(44) 



Now we substitute the expression for the kernel Eq. (|i~8|) and Eq.(|D]) into Eq. 
Using the same technique as used for Eqs.(|2"l|), (^5]) and fl2"E|), we find that there is 
only the structure B^iv )TB mil (v ) on the right hand side of Eq.(fP]). Substituting 
the explicit expressions for V\ and V 2 in Eq. Q3T)|) , and using the integration formula 
in Appendix A, we arrive at the following expression for £(1) after some tedious 
calculations 



£(i) 



Pt^Pt 



Mp 2 t)h(\Pt\) - \p 2 t C{p 2 t )h 2 (\p t 



--P 2 t D(Pt)h3(\ P t\) + ^rv 2 M\p t \)h A (\p t \) 



4vT 2 E + m D -W p 
1 , 1 

I'd 

where hi(\pt\)(i = 1,2,3,4) are given by the following equations 



(45) 



MM) = J^{^KF 1 (\pl\ qt \)[A( q 2 )-- q 2 D( q 2 )} + 



3(Q? - /i 2 ) 

[F 2 (\p t \, \q t \,fi) - F 2 (\p t \, IqtlQ^WpAiq 2 ) - ^q 2 C(q 2 ) 

-\w p q 2 b{q 2 t )] - SvTKiMH, \q t \)[A(p 2 ) - Ip 2 d(p 2 )]}, (46) 
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MN) = / ^ L,^. ^ [M\pt\,h\^)-F 2 (\p t \,\q t \,Q 1 )]A( q f) 



4tt 2 l 3(Q 2 -/i 2 
+87T«F 3 (|p t | > \q t \)\c(q 2 ) + -^l—[F A {\p t \, \q t \,») 

-^(bti, |g*i, Qi)]^fe 2 ) - STivcFiflptl, |g*|)C(p 2 )}, (47) 
Pi 



^(Iftl) = / ^{STr/cfidftUftl)^) + J,!^ J F 2 (kUg t U) 



4vr 2 — — — 3(Q?-// 



p 2 ™> 3(Q 2 - y u 2 ) 

Pi 

-87ricFi(|ft|, \qt\)[-A(p 2 t ) +p 2 t D(Pt)}}, (48) 

Mlftl) = / ^ U(Qj-U mptl lqtlfi) " Fa(N; N ' 

+87r«F 8 (lft|, k.D^C'Cg, 2 ) - ^ 2 J F 4 (M, A*) 
Pt 3 (<5i - A* ) 

M, Qi)]- 2 D{q 2 t ) - ^KF^pl \q t \)C(p 2 )}. (49) 
Pt 

The B-S scalar functions A(p 2 ), C(p 2 ) and D{p 2 ) should be normalized such that 
they satisfy Eg . ([45|) . 

V. Numerical solutions for the B-S wave function 



In this section we solve the three coupled integral equations, Eqs.(|3TD, ( p8|) and 
(|39|), numerically. The method is to discretize the integration region into n pieces 
(with n sufficiently large). In this way, the integral equations become matrix equa- 
tions and the B-S scalar functions A(pf), C(p 2 ) and D(pf) become n dimensional 
vectors. The matrix equations obtained in this way can be written in the following 
form, 

A = ZxA+ Z 2 C + Z 3 D, (50) 
R X A + R 2 C + R 3 D = 0, (51) 
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TxA + T 2 C + T 3 D = 0, (52) 
where Z i: R i: %{% = 1, 2, 3) are n x n matrices and are given by Eqs. (|3"7D, ( |3"8D and 

Substituting Eqs. (|5l"l) and (|52"D into (|50|) we obtain the eigenvalue equation for A 

#i = A, (53) 

where if is an n x n matrix 

H = Zi+Z 2 (T 3 _1 T2— -R3 1 R 2 )~ 1 (R 3 l R\— T 3 1 Ti)+Z 3 (T2 1 T 3 — Rs)^ 1 (R-2 1 Ri—T 2 x Ti) 

(54) 



The eigenvalue equation (|53D is solved by the so-called inverse iteration method 
In this way, we first construct the operator 



where A is some parameter which is chosen to be near to the eigenvalue 1 in Eq. (|53|) . 
In order to solve for the eigenvector A, we start with an arbitrary vector Y and 
operate K on Y sufficiently many times so that the eigenvector corresponding to 
the eigenvalue 1 dominates. In this way, the scalar function A is obtained. 

In our model we have several parameters, Ql, m D and E . The param- 

eter Q\ has been described in Section III, with Q\ = 3.2GeV 2 from the data of the 
electromagnetic form factor of the proton. It is noted that this value corresponds 
to the (qq') axial vector diquark (q, q' = u or d), i.e., for Eg . In the cases of Eg 
or tig- 1 this value might be somewhat different because of SU(3) flavor symmetry 
breaking. However, we do not have data to extract Q\ for Hq^ and Q^q 1 at present. 
In this work, we simply use the same value for Q\ based on the approximate 577(3) 
flavor symmetry. On the other hand, the binding energy should satisfy the following 
relation 

m w =m Q + m D + E , (56) 

Q 
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where we have omitted corrections of 0(1/i71q), since we are working in the heavy 
quark limit. Note that mr, + Eq is independent of the flavor of the heavy quark, 
because of the SU(2)f x SU(2) a symmetry. From the B-S equation solutions in the 
meson case, it has been found that the values mj = 5.02GeV and m c = 1.58GeV 



give predictions which are in good agreement with experiments [fLq| . Hence in the 
baryon case we expect 

m D + E = 0.88GW(for E^), 1.07GeF(for 1.12GeV(for (57) 

The parameter cannot be determined, although there are suggestions from 
the analysis of valence structure functions that it should be around 0.9GeV for non- 
strange diquarks [19|]- Hence we let it vary within some reasonable range. When 



we solve the eigenvalue equation Eq . (|53|) , the condition that the eigenvalue is 1 
provides a relation between a se fj and k. As discussed in Ref.|| k is related to k' (V 
is the confinement parameter in the heavy meson case and is about 0.2GeV 2 [14], III 



where k ~ Aqcd k '- Therefore, in our numerical calculations we let k vary in the 
region between 0.02GeV 3 and O.lGeV 3 . The diquark mass, mo, is chosen to vary 
from 0.9GeV to lGeV for from l.lGeV to 1.2GeV for E ( q\ and from 1.15GeV 
to 1.25GeV for vt§ . Then we obt ain the parameter for different values of mjj 
and k. The numerical results are shown in Tables 1,2 and 3 for £q\ Sq and £Iq\ 
respectively. 

With the parameters in Tables 1,2 and 3 we obtain the numerical solution for 
the B-S scalar function A as the eigenvector of Eq . (|53|) . Consequently, we get the 
numerical solutions for C and D from Eqs. (|5lD and (0). These solutions depend on 
the parameters mu and k. In Figs. 2, 3 and 4 we show the shapes of A, C and D 
for Eq\ Eq and Qq respectively. Figs. 2(a), 3(a) and 4(a) show the dependence 
on n for a typical mo, while Figs. 2(b), 3(b) and 4(b) show the dependence on mo 
for a typical k. It can be seen from these figures that for different heavy baryons the 
shape of the B-S scalar functions are rather similar. This arises from the approximate 



Table 1: Values of k and for Eg with three sets oi m D 



m D (GeV) 






0.90 






/<GeV 3 ) 


0.02 


0.04 


0.06 


0.08 


0.10 




0.5190 


0.5593 


0.5842 


0.6061 


0.6149 


m D (GeV) 






0.95 






/t(GeV 3 ) 


0.02 


0.04 


0.06 


0.08 


0.10 




0.5889 


0.6123 


0.6285 


0.6406 


0.6502 


m D (GeV) 






1.0 






/t(GeV 3 ) 


0.02 


0.04 


0.06 


0.08 


0.10 



0.6414 0.6560 0.6669 0.6757 0.6828 



Table 2: Values of k and for 2g with three sets oi m D 



m D (GeV) 






1.10 






/t(GeV 3 ) 


0.02 


0.04 


0.06 


0.08 


0.10 




0.5047 


0.5402 


0.5643 


0.5826 


0.5974 


m D (GeV) 






1.15 






/«(GeV 3 ) 


0.02 


0.04 


0.06 


0.08 


0.10 




0.5785 


0.5995 


0.6155 


0.6283 


0.6391 


m D (GeV) 






1.20 






/t(GeV 3 ) 


0.02 


0.04 


0.06 


0.08 


0.10 



0.6341 0.6478 0.6588 0.6682 0.6763 



Table 3: Values of k and for f2g with three sets of m D 



m D (GeV) 






1.15 






/<GeV 3 ) 


0.02 


0.04 


0.06 


0.08 


0.10 




0.4975 


0.5325 


0.5565 


0.5748 


0.5897 


m D (GeV) 1.20 


/t(GeV 3 ) 


0.02 


0.04 


0.06 


0.08 


0.10 




0.5729 


0.5935 


0.6093 


0.6222 


0.6331 


m D (GeV) 






1.25 






«(GeV 3 ) 


0.02 


0.04 


0.06 


0.08 


0.10 




0.6296 


0.6430 


0.6539 


0.6633 


0.6714 
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SU (3) flavor symmetry and is to be expected. All the scalar functions decrease to 

zero when \p t \ is larger than about 1.5GeV, because of the confinement interaction. 

(*) 

Furthermore, since we are discussing the ground states ujq , there are no nodes in 
the functions. 

VI. Application to the nonleptonic decays — ► fiW P (V) 

In this section we will apply the numerical solutions of the B-S equation to the 
nonleptonic decays Qb ^c*"* an d a pseudoscalar or vector meson. In fact, 
and Ej^ decay strongly and their weak decays are hard to observe. However, Qj, 
decays only weakly. We will first calculate the Isgur-Wise functions and C( a; ) 
for — > in Eq.(fTO|) and then apply them to the nonleptonic weak decays of 

n b . 

A. Isgur-Wise functions for — > fiW 

The Isgur-Wise functions £(o>) and CO^) are related to the overlap integrals of 
the B-S wave functions of the initial and final (fiM) states. The concrete 
expression for them can be obtained by comparing the structure B^,(v')TB miM (v) 
and v ■ B m i{v')Tv' ■ B m (v) on both sides of Eq. (|40"D . Similarly to Eq. (|44|) , we have 
the following equation after substituting Eq.(^) into Eq.(pE2"|) and using Eq.( 



d 4 q 



^ 2n ^G^(P,p,q)xp m (p). (58 
Substituting Eq.(|17]) and Eq.flU) into Eq.® and using Eqs.(g§, (g5|) and (||) 



we find that on the right hand side of Eq. (|58|) there are the following structures: 
p\ ■ B m i(v')Tv' ■ B m (v), v ■ B m i(v')Tp t ■ B m (v), p' t ■ B ml (v')Tp t ■ B m (v), p t ■ B m ,(v')Tp t ■ 
B m (v), and p' t ■ B m i{y')Tp' t ■ B m {v). However, all of them can be expressed in terms 
of B^,{v')TB mii (y) and v ■ B m i{v')Yv' ■ B m (v), after the integration over p, on the 
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grounds of Lorentz invariance. Take p' t ■ B m i(v')Tp t ■ B m/I (v) as an example. In 
general, the integral / jk^sp'tPtfi where / is some Lorentz scalar function, can be 
expressed in terms of g^, v^v v , v'^v'", v'^v v and v^v' v . However, only the 
and v^v' u terms contribute when contracted with B m in{y')TB mv (y), leading to the 
structures B 1 ^, {v')Y B mil {v) and v ■ B m /(v')Tv' ■ B m (v), respectively. The coefficients 
of these two terms can be obtained directly. In this way, we have the following 
replacement rule: 

p'fB m ,(v') r 
v ■ B ml (v') r 
p'fB m ,(v') r 

+ 

p t -B m >(v') r 
+ 

p' t -B m r(v') r 
+ 

Since in the weak transition the diquark acts as a spectator its momentum in 
the initial and final baryons should be the same, P2 = p' 2 - Then we can show that 

p' = p + m D {y ' — v), (60) 

where we have omitted the 0{l/mo) corrections. From Eq. ( |B"0"D we have the 
following relations straightforwardly 

p[ = piuj - ptVoo 2 - lcos6>, 

p t 2 = p\ + p\ {u 2 — l)cos 2 # + p 2 (uj 2 — 1) — 2pip t uVuj 2 — lcOS^, 
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v' ■ B m (v) -> 

1 — uj z 

1 



v ■ p' t v ■ B m f(v')Tv' ■ B m (v) 



p t ■ B m (v) -> - v' ■ p t v ■ B m >(v')Tv' ■ B m (v), 
1 — uo z 



p t ■ B m (v) 



-. — - -p t ■ p. H — t— —v ■ p.v • Pt 

2(uo 2 -\y yt 2(uj 2 -!) 2 n yt 



1 , UJ 

2 Pt ' Pt ~ 2(uj 2 -1) 

2 j 2 



v ■ p' t v' ■ pt 



BUv)TB m ^v 



p t ■ B m (v) 



1 2 1 

~2 Pt 2(uj 2 - V 



v ■ B m ,{y')Yv' ■ B m (v) 
Bt,(v')TB mfl (v) 
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3uj 



2(uj 2 - 1 

Pt ■ B m( V ) 
UJ 

2(uj 2 - 1 



2(uj 2 - 1) 



;(v' -Ptf 



v ■ B m ,{v')Tv' ■ B m (v) 



i /2 | 1 

2 Pt 2(uj 2 - 1) 



/\2 



■ Pt) 



B^(v')VB mil (v 



3uj 
2(uj 2 - 1) 



(v-p't) 2 



vB m ,{v')Tv' -B m {v). (59) 



v ' ■ pt = —ptV — lcos^, 
v ■ p' t = (1 — uj 2 )pi + up t Voo 2 — lcos#, 

Pt ■ p't = —Pt — Pti.u 2 — ^)cos 2 9 + uj\Juj 2 — lpip t cos6, (61) 

where 9 is the angle between p t and v' t . 

With the aid of the relations between A(p), C(p), D{p) andv4(p 2 ), C(p 2 ), D{p 2 ) 
[Eqs.(|5^)-([5^)1 and using Eqs.(^), fl6"T|) and the integration formulae in Appendix A 
we have the explicit expressions for £,{uj) and (,(oS) after integrating the pi component 
by selecting the proper contour, 

£H = l^r^m ± == 

J 4n 2 Jo 2W P (E + m D -W p )(E + m D - uW p - p t Vuj 2 - lcos#) 

{-2Wy£ + m D - Wgi^fo 2 , cosfl)i(p 2 ) - - cos 2 fl)MM) - /i 3 (b t |)] 



F A (pl cos6) - -ujp 2 (l - cos 2 6)F c {pl cos9)h 2 (\p t \) + -p^u 2 - 1(1 - cos 2 £) 
cosflMH) " MbtDJi^CP?, cos #) + ^i 1 - cos 2 6)(E + m D - W p )p 2 A(p 2 ) 

F D (p 2 , cos9) - l -p\{\ - cos 2 #)[(^ 2 - 1)W P + LuptVuj 2 - lcos6]h 2 (\p t \) 
3 



F D (p 2 , cos6) + -(1 - cos 2 fl)[p 2 + p t 2 (cj 2 - l)cos 2 # + ojVoj 2 - lW p p t ] 
Fniplcose^dptD-hsdptl)]}, (62) 



and 



cm = r ^ - i— 

J 4vr 2 7o 2 Wp(£ + m D - W p ) (£ + m D - uW p - p^uj 2 - lcosfl) 

{1 3cj 
-== Pt cos0F A (p^cos0)fc 2 (|p t |) + 2 ^ (1 - 3cos 2 #) 
— 1 4 (or — 1 ) 

M\Pt\) -h(\p t \)]F A (p 2 t , cose) + W p + ^— Pt cos9 F c (p 2 ,cos9) 

yoj z — 1 

2W P (E + m D - W p )A(p 2 ) + J_ - [iup 2 - 3^cos 2 % 2 

3 

-2W p p t \/oj 2 - lcos6]F c (p 2 ,cos6)h 2 (\pt\) - —j= -- [copt - Sucos 2 9p t 

4\/uj z — 1 

-2^ v / ^Tcos0]cos0[/i 1 (|^|) - h 3 (\p t \)}F c (p 2 t ,cos9) + \-uW 2 
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[-2W P {E + m D - W p ))F D {p 2 t , cos9)A(p 2 ) + - X^— Wt ~ 3^cos 2 

zyw — 1 



-2W p p t VuJ 2 - 1cos^](v / cj 2 - 1W P + ujp t cos6)F D (p 2 , cos6)h 2 {\p t \) 
3 



~A(uu 2 - 1) 



[uip t — 3ucos 2 9p t — 1W v \Joj 2 — lcos#][p t +p i (c<j 2 — l^cos 2 ^ 



+^v / ^TW^cos^]F i ,(p 2 ,cos^)[/ il (|^|)-/i3(bt|)]}, (63) 



where = 1, 2, 3, 4) are given in Eqs. (|46l) -([49|) and Fa{p 2 , cos8), F c (p 2 , cos8 

and Fc(p 2 ,cos9) have the following expressions 



F A (p 2 ,cos£) = 8vr K F 1 (|KUg,|)[i(g 2 )--g 2 D(g 2 )-i(^ 2 )] + 



3(Q 2 -/i 2 ) 



G% 2 ) - ^(% 2 )) - \q 2 M q 2 )} + ^ 3(Q 1 f 6 ^ 2) [-F 4 (b;i i 
\qt\,Qi)][C(q 2 ) + 2(cuW p + p t V^lcos9)D(q 2 )] 

1 



+^8™F 5 (|p;Uft|)D(g?)J, (64) 
F c (p 2 ,cos£) = ±J^{itoid?iMM)[{u>W p +^ 

167T/5 



+p'tD(p' 2 )) + {{uW p + p t Vu^lcos9f - m 2 D )C(p' 2 )} 



3(Q? ~ V 2 ) 



[F 2 (| P ;|, |g t |,/i) - F 2 (| P ;|, |g t |,gi)][(wWp + ^v / ^Icosfl) 2 + m 2 D ] 

Mtf) ~ -^StckFMI, \q t \)[(cuW p + p t V^lcos6) 2 - m 2 D ]C(q 2 ) 
Pt 

-— 16n(3 <-F 4 (\p[\Mri + FM\Mt\,QiM"W p 



v? 3(g 2 - n 

+ Pt V^ 2 ^lcos9)C(q 2 ) + {{uoW p + p tV / ^ 2 ~ 3 Icos#) 2 + m 2 D )L)(q 2 t )} 
-^toKFMl \q t \)(ujW p + p t V^T C ose)D(q 2 )\ , (65) 



F D (p 2 ,cos£) = J ^ 18^(1^1,1^1) 



A^) + 3?^(% 2 ) - MP?) 
z Pt 



+(uW p + ^v^^lcose)^^ 2 ) + (pf + m 2 D )D(p' t 



r2\ 



167T/3 



3(0? - ^ 2 ) 



[F 2 (K|, \q t \,fi) - F 2 (\p' t \, \qt\M] (uW p + p t V^lcos6) (A(q 2 ) 



2 \ 2 

m o Jn/J\ I i m D„2ni„2\ 



^ D{qf) + ^gg C(q I 

Pt J m 
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-1^87TKF 3 (\p' t \,\q t \)(uW p 
Pt 



+p t Vuj 2 - lcos0)C% 2 ) - -4 Q /n2 67r/9 2 \ l~ F ^P'^ 1*1'^) 



+^(1^1, |g t |,QO] 



^|±M(C7(g 4 2 ) + 2(u,W p + p t V^lco S 0) 



L>(g 2 )) - (cuiy p + p iV / ^Icos^) J D(g 2 ) 



3m 2 ) + 2pf 



toKFM,\q t \)[D(q*)}. (66) 

In Section V we obtained the numerical results for A(pf), C(pf), and D(pf). 
Substituting these results into Eqs.(|62|) and (|63|) we have the numerical solutions 
for £(u) and ((uj) depending on the parameters in our model. For fij,*^ — > fiM, 
we show the Isgur-Wise functions with typical value = 1.20GeV in Fig.5(a) 
(k = 0.02GeV 3 ) and Fig. 5(b) (k = O.lOGeV 3 ) respectively. The dependence of the 
Isgur-Wise functions on mo are shown in Fig. 5(c) (mjj = 1.15GeV) and Fig.5(d) 
(rriD = 1.25GeV) for /t = 0.06GeV 3 . It can be seen from these plots that £(^) and 
C(to>) have opposite signs and £(cj) changes more rapidly than ((u) as u increases. 

It is interesting to study the relation between £(u;) and ((uj). Based on the 
picture that in the large N c limit heavy baryons are viewed as the bound states of 
chiral solitons and heavy mesons [PI, Chow has shown that £(co>) and ((u>) obey the 
following relation PT[] 

e(w) = -(l+w)C(a;). (67) 

The deviation from this relation is caused by 1/N C corrections. From Figs.5(a)-(d) 
we can see that, in the range of the parameters in our model, Eq. floTf) is generally 
satisfied. For some sets of parameters this relation holds well. 

B. Nonleptonic decays Q 6 -> fiW P (V) 

In this section we will discuss the Cabbibo-allowed two body nonleptonic decays 
of ilb ^c*"* P (V) (P an d V stand for pseudoscalar and vector mesons respec- 
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tively). The Hamiltonian describing such decays reads 

H cS = 9± VcbV * D ( ai0l + a 2 2 ), (68) 

with 0\ = (DU)(cb) and 2 = (cU)(Db), where U and D are the fields for light 
quarks involved in the decay, and (51^2) = <Zi7,u(l — 75)92 is understood. The pa- 
rameters ai and a 2 are treated as free parameters since they involve hadronization 
effects. Since Q b decays are energetic, the factorization assumption is applied so 
that one of the currents in the Hamiltonian QB"8] ) is factorized out and generates a 
meson 



|22| , [23| . Thus the decay amplitude of the two body nonleptonic decay be- 
comes the product of two matrix elements, one is related to the decay constant of 
the factorized meson (P or V) and the other is the weak transition matrix element 



between Qh and &i*> 



M^(n b - fiWp(V)) = ^K^ D a 1 (P(\/)|A M (^)|0)(fiWK)|Jl^M). (69) 

Here J M is the weak current (cb) and (01^(1/^) \P(V)) are related to the decay 
constants of the pseudoscalar or vector mesons by 

(0K|P> = ifpq^, 

W^V) = fvmve^ (70) 

where g M is the momentum of the emitted meson (from the W-boson), e M is the po- 
larization vector of the vector meson, and the normalization for the decay constants 
is chosen so that = 132MeV. It is noted that in the two body nonleptonic weak 
decays tt b -> fiWp(V) there is no contribution from the a 2 term, since such a term 
corresponds to the transition of Q b to a light baryon instead of . On the other 
hand, the general form for the amplitudes of Q b QyP(V) are 

M(Q b ^Q c P) = iu f (v')(A + B l5 ) Ui (v), 

M(tt b ^tt c V) = u f (v')e* fl [A 1 ^ 5 + A2Pf^5 + Bi^ + B 2 p f Aui(v), 
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M(n b ^tt* c P) = i qtl u>}(v')(C + D l5 ) Ui (v), 

M(Q b ^Q* c V) = u u f (v')e*' l [9u,(C 1 + D ll5 )+p tu ^(C 2 + D 2l5 ) 

+PivPfn{C 3 + D 3 j 5 ))ui(v), (71) 

where is the dirac spinor of il b , uf' is the dirac (Rarita-Schwinger) spinor of 
and is the momentum of fl b (f^*)). 

From Eqs.(|69|)- ([7T|) and using Eq . (f40|) we find 

a = ^k^^i/p - %) [(^ + 2)e(w) + (w 2 - i)ch], 

B = ^V cb V{} D a 1 f P ^(m l + m f )[(3u-2)auj) + 3(u 2 -l)auj)}, 
A, = B 1 = ^V^ D a 1 f v m v ^[Lo^) + (u J 2 -l)C^)}, 
A 2 = ^V^aJvmvU— -— ) [£(«) + (« + l)C(u)], 

B 2 = -^V A V^ D a 1 fvm v U— + —)[^u) + (u-l)acj)], 
V 2 o \rrii wif j 

c = -^K^Vi/p-^(i + ^-)[eM + (^-i)CM], 

G F 1 
D = --^VcbV^aJp-^ 

G 2 
Ci = Di = -^ ch Vu D a x f v m v -j=£{uj), 

C 2 = -^V cb V* Dai f v m v ^= — [f( w ) + ( w +l)CM], 
V 2 v 3 rrii 

D 2 = ^Lv^aJvmv^ — [£{u) + {u - l)((u)], 

V 2 V <J m i 

C 3 = D^I^V^^CM, (72, 

where (mj) is the mass of £l b (f^). 

With Eqs. (|7TD and (|72|) we can calculate the decay widths and polarization pa- 
rameters for Cl b — > P (V). The kinematic formulae which have been derived 



^ m f 

rrii 



)e(w)+(w-l-(a; + 3)^)c(a;) 



using both partial wave and helicity methods can be found in references [24, 25 



These two methods are equivalent. For instance, in the helicity method[25|, the 
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decay width is expressed in terms of the helicity amplitudes, 



r " 16^? ^ l^poo^l 2 . (73) 



* A,;, A 



i,A/ 



where p c is the cm. momentum of the decay products and the helicity amplitudes 
are defined as 

hx fMvy ,x t = {n£\\ f ),P(y)(\pM)\H*MXi)) (A/ - X P(V ) = Ai). (74) 
The "up-down" asymmetry is given by 

a _ Y,\ f {\h\ f .\ p(v) -l/2\ 2 ~ I^Ay.Aptv-);-!^! 2 ) ^ 



The relations between the helicity amplitudes and the amplitudes given in Eq.([7T|), 
which we will not list here, can be found in |25|, EfJ. Then from Eqs.([T2)-([75D, we 
obtain the numerical results for the decay widths and asymmetry parameters. In 
Table 4 we list the results for = 1.20GeV. The numbers without (with) brackets 
correspond to « = 0.02GeV 3 (k = O.lOGeV 3 ). The results for k = 0.06GeV 3 in 
the range mo = 1.15GeV (without brackets) and mp = 1.25GeV (with brackets) 
are shown in Table 5. In the calculations we have taken m^ b = 6.14GeV and the 
following decay constants 



U = 132MeV, f Ds = 241MeV27, f p = 216MeV, f Ds = /„.. 



It can be seen from Tables 4 and 5 that the predictions for the decay widths show 
a strong dependence on the parameters k and mo in our model. The experimental 
data in the future will be used to fix these parameters and test our model. However, 
the dependence of the up-down asymmetries on these parameters is slight. 

The decay widths and asymmetry parameters have also been calculated in the 
nonrelativistic quark model approach |2B]], where the form factors are calculated at 
the zero-recoil point and then extrapolated to other uj values under the assumption 
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Table 4: Predictions for decay widths and asymmetry parameters for — > 
fiW p (y) for m D = 1.20GeV. 



Process r(10 10 s x ) a 

fy -> QP c u- 0.052a? (0.154a?) -0.67 (-0.70) 

fy -> fi°L>7 0.261a? (0.592a?) -0.56 (-0.58) 

-> n»p~ 0.073a? (0.207a?) -0.68 (-0.71 J 

fy -> fi^DT 0T15oi (0.245a?) -0.73 (-0.74) 

fy ->■ fi;°7r- 0.046a? (0.133a?) -0.61 (-0.58) 

fy -> nfD- 0.165a? (0.370a?) -0.54 (-0.52) 

IV ->■ QJV 0.134a? (0.354a?) 0.59 (0.59) 

fy -> 0.462a? (0.960a?) 0.31 (0.31) 



Table 5: Predictions for decay widths and asymmetry parameters for Q b 
fiW p (y) for k = 0.06GeV 3 . 

Process r(10 lu s~ 1 ) a 

fy -> QP c u- 0.075a? (0.145a?) -0.64 (-0.72) 

fy -> flgDj 0.358a? (0.562a?) -0.54 (-0.59) 

iy -> figp- 0.102a? (0.150a?) -0.65 (-0.70) 

iy ->■ fi°Zy~ 0.149a? (0.232a?) -0.71 (-0.75) 

fy ->■ 0.067a? (0.123a?) -0.64 (-0.56) 

IV ->■ figgj 0.227a? (0.345a?) -0.56 (-0.50) 

n b -> fi*V 0.200a? (0.314a?) 0.59 (0.59) 

iy -> fi;°£T 0.616a? (0.888a?) 0.31 (0.31) 



of a dipole behaviour. By comparing the predictions in the B-S and quark models 
we find that the decay widths in our model are smaller than those in the quark 
model. For the asymmetry parameters, the difference is even larger. Except for the 
processes fy — ^c* 71 "^ an d fi&" ~^ ^T^7 m Tables 4 and 5, even the signs of a in 
these two models are opposite. 

VII. Summary and discussion 

Since in the heavy quark limit the light degrees of freedom in a heavy baryon 
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have good spin and isospin quantum numbers and since the internal structure is 

blind to the flavor and spin direction of the heavy quark, we assume that a heavy 

(*) 

baryon, ojq , is composed of a heavy quark and a light axial vector diquark. Based 
on this picture, we establish the B-S equation for the heavy baryon uJq\ We discuss 
the form of the B-S wave function and find that in the heavy quark limit there 
are three B-S scalar functions to describe the dynamics inside a heavy baryon Uq\ 
This is consistent with our physical picture. In order to solve the B-S equation, 
we assume a kernel containing a scalar confinement term and a one-gluon-exchange 
term, as in the Aq case. In the heavy quark limit, the heavy quark is almost 
on mass-shell inside a heavy baryon and it is appropriate to apply the covariant 
instantaneous approximation in the kernel. Then we derive explicitly three coupled 
integral equations for the three B-S scalar functions A, C and D. These equations 
are solved numerically and we give the model predictions for these functions. The 
results appear reasonable. It is shown that the shapes of these functions are similar 
for Eq\ E^q and Qq\ with differences arising from SU(3) flavor symmetry breaking 
effects. 

Although the B-S equation is formally the exact equation to describe the bound 
state, there is much difficulty in applying it to the real physical state. The most 
difficult point is that we must take a phenomenologically inspired form for the kernel. 
Furthermore, we have used the quark and diquark propagators with their free form, 
which leads to some uncertainties. In our approach, there are several parameters 
such as k, m D and a;( eff \ subject to the condition that the observed masses of Uq* 
be reproduced. In our numerical solutions we let these parameters vary in some 
reasonable range. Another parameter is Q\, which arises from the internal structure 
of diquark. Its value for the (qq') diquark (q, q' = u or d) is extracted from the data 
of the electromagnetic form factor of the proton. When there is a strange quark in 
the diquark, we do not have a means to determine its exact value at present. In 
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the future, the experimental data for uJq should help to fix the parameters in our 
model. 

As phenomenological applications, this model has been used to calculate the the 
Isgur-Wise functions and ((uj) for — > and consequently, has provided 
theoretical predictions for the Cabbibo-allowed two body nonleptonic decay rates 
and up-down asymmetries for the physical processes fi 6 — > P (V). It has been 
shown that the relation between £(a>) and ((u) in our model is generally consis- 
tent with that in the soliton model in the large N c limit. We have also compared 
our model predictions with those in the nonrelativistic quark model. Our model 
yields decay widths which are much smaller and for the asymmetry parameters the 
difference is even bigger. All these predictions will be tested in future experiments. 
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Appendix A. Integration formulae 

In this appendix we give the formulae which are used to reduce the three di- 
mensional integration to the one dimensional integration. In the following formulae 
4>{q 2 ) is some arbitrary function of qf. The relevant results needed are: 

with 

2 

F i(\Pt\, M) = 7-2- — — t — ^ — 7-3-2 ■ (77) 
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(2tt)3 (p t - qt )2 + 52 

with 



' l?,U) = «R ln (N-k,l) 2 + ^' (79) 



d 3 g t Pt • gt0(g f 2 ) 



47T 2 



#£)*s(|ft|,|ft|), 



JO) 



with 



*3(|ft|,|ft|) 



4btlkt| 



In 



(|Pt|-|&l) 2 + /^ , MPt\\Qt\(Pt + Qt 



+ 



(\Pt\ + \<h\) 2 + H 2 (p 2 t+q 2 t+V 2 ) 2 -±vWt 



*1) 



/ 4 = 



d 3 g t (pt ■ qtfHq 2 ) 

(27T) 3 (p t - % ) 2 + 5 2 



-/ 



4tt 2 



0(? t 2 )^(bt|,|? t |,5), 



(82) 



with 



^4(1^1,1*1,5) = 



p 2 + q 2 + ^ 2 



x , p 2 + gf + a 2 ln (N 



|g f |) 2 + ^ 2 



4|ft||ft| (M + |g t |) 2 + 5 2 



(83) 



d 3 * (pt-ft)W) 

(27T)3[(p t - ?t ) 2 +^ 2 ] 2 



/ 



g 2 dg f 

4tt 2 



0(g 2 )F 5 (|p t |,M), 



(84) 



with 



^6(|ft|,|ft|) = 2 



1 + 



p\ 



+ q 2 t + »\ n (\Pt\ - \qt\) 2 + /x 



(P? + 9t + 



2|ft||ft| (\Pt\ + \qt\) 2 + V 2 ' (p? + q 2 t + V 2 ) 2 - 4p?g? 



(85) 
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Figure Captions 



Fig.l The vertex of two axial vector diquarks and a gluon. 

Fig. 2 The B-S scalar wave functions for E^. The units are: GeV~ 2 for A(p 2 ), 
GeV~ 3 for C(p 2 ) and GeV" 4 for D(p 2 ). (a) For m D = 0.95GeV, we show the de- 
pendence on \p t \ for two values of k. In the upper plane, the upper (lower) solid 
line is for A(p 2 ) (C(p 2 )) with k = 0.02GeV 3 , the upper (lower) dotted line is for 
A{p 2 ) {C(p 2 )) with k = O.lOGeV 3 . In the lower plane, the solid line is for D{p 2 ) 
with k = 0.02GeV 3 and the dotted line is for D(p 2 ) with k = O.lOGeV 3 . (b) For 
k = 0.06GeV 3 , we show the dependence on \p t \ for two values of mp. In the upper 
plane, the upper (lower) solid line is for A{p 2 ) {C{p 2 )) with m D = 0.90GeV, the 
upper (lower) dotted line is for A{p 2 ) (C(p 2 )) with m D = l.OOGeV. In the lower 
plane, the solid line is for D(p 2 t ) with m D = 0.90GeV and the dotted line is for 
D(p 2 ) with m D = l.OOGeV. 

Fig. 3 The B-S scalar wave functions for Sq\ The units are: GeV -2 for A(p 2 ), 
GeV -3 for C(p 2 ) and GeV -4 for D{p 2 ). (a) For m D = 1.15GeV, we show the de- 
pendence on \p t \ for two values of k. In the upper plane, the upper (lower) solid 
line is for A(p 2 t ) (C(p 2 )) with k = 0.02GeV 3 , the upper (lower) dotted line is for 
A(p 2 ) {C(p 2 )) with k = O.lOGeV 3 . In the lower plane, the solid line is for D(p 2 ) 
with k = 0.02GeV 3 and the dotted line is for D(p 2 ) with k = O.lOGeV 3 . (b) For 
k = 0.06GeV 3 , we show the dependence on \p t \ for two values of m D . In the upper 
plane, the upper (lower) solid line is for A(p 2 ) (C(p 2 )) with m D = l.lOGeV, the 
upper (lower) dotted line is for A(p 2 ) (C(p 2 )) with m D = 1.20GeV. In the lower 
plane, the solid line is for D(p 2 ) with = l.lOGeV and the dotted line is for 
D(p 2 ) with m D = 1.20GeV. 

Fig. 4 The B-S scalar wave functions for f^. The units are: GeV~ 2 for A(p 2 ), 
GeV" 3 for C(p 2 ) and GeV" 4 for D(p 2 ). (a) For m D = 1.20GeV, we show the 
dependence on \p t \ for two values of k. In the upper plane, the upper (lower) solid 
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line is for (C(pj)) with k = 0.02GeV 3 , the upper (lower) dotted line is for 

A(pf) (C(pf)) with k = O.lOGeV 3 . In the lower plane, the solid line is for D(pf) 
with k = 0.02GeV 3 and the dotted line is for D(pf) with « = O.lOGeV 3 . (b) For 
k = 0.06GeV 3 , we show the dependence on \p t \ for two values of mo- In the upper 
plane, the upper (lower) solid line is for A{pf) (C(pj)) with m D = 1.15GeV, the 
upper (lower) dotted line is for A(pf) (C(pf)) with m D = 1.25GeV. In the lower 
plane, the solid line is for D(pf) with m D = 1.15GeV and the dotted line is for 
D(p 2 t ) with m D = 1.25GeV. 

Fig.5(a)-(d) Numerical solutions for and ({uj) for — > The upper 
solid line is ^(o;) and the lower solid line is C(^)- The dotted line is + 1). 

The parameters are mo = 1.20GeV and k = 0.02GeV 3 in (a), m D = 1.20GeV and 
k = O.lOGeV 3 in (b), « = 0.06GeV 3 and m D = 1.15GeV in (c), « = 0.06GeV 3 and 
m D = 1.25GeV in (d). 
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Fig.3(b) 




Fig.4(b) 
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Fig. 5(a) 




1 1.2 1.4 1.6 1.8 2 2.2 



Fig.5(b) 
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Fig. 5(c) 
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Fig.5(d) 



